Effect of sigma meson on the 1)1(2430) — )■ Dtctt decay 
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We study the effect of signia meson on the Di (2430) — >■ Dnn decay by constructing an effective 
Lagrangian preserving the chiral symmetry and the heavy quark symmetry. The sigma meson 
is included through a linear sigma model, in which both the qq and qqqq states are incorporated 
respecting their different U{1)a transformation properties. We first fit the sigma meson mass and a- 
TT-TT coupling constant to the / = 0, S'-wave tt-tt scattering data. Then, we show how the differential 
decay width dr(_Di — > D{nn)i^o,L^o)/dmT„T depends on the quark structure of the sigma meson. 
We find that our study, combing with the future data, can give a clue to understand the sigma 
meson structure. 

PACS numbers: 14.40.Rt, ll.30.Rd, 12.39.Hg 



I. INTRODUCTION 

The lightest scalar meson "sigma" is an interesting ob- 
ject which may give a clue to understand some funda- 
mental problems of QCD such as the chiral symmetry 
structure, the origin of mass and so on. The mass spec- 
trum of the light scalar meson nonet including the sigma 
meson disfavors the qq picture but prefers the qqqq inter- 
pretation (see, e.g. Ref. [ll, 0] and references therein). 

If in the nature there are both qq and qqqq scalar states, 
they might mix to give the physical scalar mesons (some- 
times the glueball component is included). In the liter- 
ature, this superposition has been discussed widely Q, 
but the structure of the light scalar mesons remains an 
open question and deserves further investigation. 

To investigate the quark contents of the light scalar 
mesons, we should find some quantities that can distin- 
guish different components, two-quark component, four- 
quark component and glueball. In Ref. [J] it was pointed 
that, in a linear sigma model expressing the spontaneous 
chiral symmetry breaking, the qq and qqqq bound states 
have different charges of U{1)a symmetry therefore one 
can use this ^'^(l)^ transformation property to discrim- 
inate the qq and qqqq components of the scalar mesons. 
The U{1)a symmetry is explicitly broken by anomaly in 
QCD which causes a mixing between the qq states and 
the qqqq states p-Ql- This contribution is included in 
the model Lagrangian in such a way that the anomaly 
matching condition is satisfied. Then, unlike the other 
terms which are invariant under the U{1)a transforma- 
tion, the U{1)a violating terms are constrained by the 
anomaly. Another source of the mixing between qq and 
qqqq comes from the existence of both the qq and qqqq 
condensation @, Q- 

In this paper, we devote ourselves to study the sigma 
meson structure in the heavy-light meson decay, specif- 
ically, the I?i(2430) — >■ D{Tnr)i=o,L=o decay, based 
on the chiral partner structure between {D*,D) and 
(1)1(2430), 1)^(2400)) 0,0 (in the following, we simply 
denote i:»i(2430) and D^{2400) as L>i and D^, respec- 
tively). An interesting property of this process is that 
there is only one light quark in the heavy-light meson so 



that the axial transformation property is well determined 
and the heavy-light meson couples only to the two-quark 
component of the scalar meson. 

The paper is organized as follows: In sec. |TT1 we in- 
troduce an extended linear sigma model for three flavor 
QCD. Section Hill is devoted to study the tt-tt scattering: 
We determine the cr- tt-tt coupling and sigma meson mass 
by fitting them to the data for the / = 0, S'-wave channel. 
In sec. IIVI we construct an effective Lagrangian for the 
interaction among light mesons and heavy-light mesons. 
In sec. |Vl after determining the relevant parameters in 
the heavy-light meson sector, we show how the differen- 
tial decay width dT{Di — > _D(7r7r)/=o,L=o)/'^"^7r7r depend 
on the quark structure of the sigma meson. We give a 
summary and discussions in sec. I VII In appendices we 
show some details of derivations of several formulas. 



II. LINEAR SIGMA MODEL WITH 
TWO-QUARK AND FOUR-QUARK STATES 

In this section, we introduce a linear sigma model for 
three flavor QCD in the low-energy region by including a 
3x3 chiral nonet field M representing the qq states and 
a 3 X 3 chiral nonet field M' standing for the qqqq states. 
These two nonets have the same chiral SU{3)l x SU{3)ji 
transformation property: 



M -^ gLMgl , M' -^ QlM'q^ , 



(1) 



where gi^ G SU{S)l,r- On the other hand, as pointed 
in Ref. [J, Q, they have the following different U{1)a 
transformation properties: 



M 



e^'^M , 



M' 



-> e 



'M', 



(2) 



with a as the phase factor of the axial transformation. 
We decompose M and M' as 



M = S + i-P , M' = S' + i$' , 



(3) 



where S and S" are the scalar meson matrices and $ and 
$' are the pseudoscalar meson matrices. 



In this paper we adopt the foUowing extended hnear 
sigma model Jj]: 



^lif 



ht 



-Vo{M,M')-VsB , 



(4) 



where the first two terms are the kinetic terms for the qq 
and qqqq states, Vq{M, M') is the potential term invari- 
ant under the chiral SU{3)l x SU{3)ii transformation, 
and VsB stands for the explicit chiral symmetry break- 
ing terms due to the current quark masses. It should be 
noted that Vo{M,M') is decomposed as 



Fo(M, A/') = VInv(M, M') + K,(M, M'), 



(5) 



where Vinv is invariant under the U{1)a transformation 
while Vri violates the U{1)a symmetry explicitly due to 
the anomaly. We do not specify the form of Vq in our 
present analysis but assume that this model allows for 
spontaneous chiral symmetry breaking by a consistent 
choice of the parameters in Vq. Since the U{1)a symme- 
try is explicitly broken by the anomaly, the form of the 
U{1)a violating term Vi, is constrained by the anomaly 
matching condition. 

The U{1)a symmetry breaking by anomaly causes a 
mixing between the qq states and the qqqq states, as 
pointed in Refs. p-Q- In. addition, the spontaneous chi- 
ral symmetry breaking also generates a mixing between 
the qq states and qqqq states ISJ, (7|. As a result, four 
physical iso-singlet scalar mesons are given as the mixing 
states of two qq states and two qqqq states through the 
mixing matrix Uf as 

(f//)la {UfU (Uf),, {UfU 
{Uf)2a {UfU (Ufhc {Uf)2d 

{Ufha {UfU {UfU {UfU 
{UfU {UfU {UfU (UfU 





values (VEVs) of the qq and qqqq scalar fields. In the 
chiral limit, we have 



where /pi, . . ., /p4 are the physical scalar mesons with 
mass ordering mf-^ < ruf^ < rn/g < rnf^. In the follow- 
ing, consistently with PDG notation, we use the notation 
a for the lightest scalar meson /pi. fa and fb are the qq 
states with fa = {Si + S'D/v^ and /& = S^ and fc and 
fd are the qqqq states with fc = {{S')\ + (S")|)/\/2 and 

ft = {S')l 

Similarly, for the two iso-triplet pseudoscalar mesons, 
we have 



Ftt cos ( 



2W2 



Ftt sin t 



-2V4 , 



(8) 



where F^ = 130.41 MeV [2] stands for the decay constant 
of 7r(140) and V2 and V4 denote the VEVs of the qq and 
the qqqq scalar fields, respectively. 

Using the method shown in Ref. [5|, we obtain the 
relation among the scalar meson- tt-tt coupling constant 
gf-rnn the mixing matrices and the scalar mass as: (for a 
derivation, see Appendix E]) 



5/,T 



F 



^{Uf)ja-sm9^{Uf)jc]ml, (9) 



which is similar to the tr- tt-tt coupling of Ref. W, except 
the mixing angle included and the chiral limit taken. We 
also obtain the relation for the four-pion coupling con- 
stant gTr-KTr-K as (see Appendix R|) 



g-K-n 



6 



4 



TP2 Z^ 

"" 3 = 1 



'cos0^{Uf)ja-sm6-^{Uf)jc] mj. 



(10) 



Making use of the relations (j9|) and ()10|) together with 
the orthonormal conditions. 



{uu=y^{u,ry^i, 



4 

E 
4 

3 = i 



4 



J = l 



U.aU.c^O: 



(11) 



we obtain the following sum rules: 



(6) 



E 



9fj7^7V 



J=l J] 



2 \2 



7 = 1 ^ fi' 



^g-mnnr j 



2 

F2 



(12) 
(13) 



We should note that, as will be shown in the next section, 
these sum rules guarantee that the tt-tt scattering am- 
plitude satisfies the low-energy theorem of the Nambu- 
Goldstone (NG) bosons. 



cos f 7r — sm f 
sin Ott cos ( 



(7) 



where TTp and tt' in the left-hand side are the physical 
states while tt and tt' in the right-hand side denote the 
qq states and qqqq states, respectively. In the present 
analysis we identify TTp as tt(140). As was shown in 
Refs. 0, [3, the above pseudoscalar mixing angle Ott re- 
lates to the pion decay constant and vacuum expectation 



III. SIGMA MESON IN THE tt-tt SCATTERING 

In this section, we determine the sigma meson mass 
Too- and cr-TT-TT coupling constant ga-Tm using the isospin 
zero (/ = 0), S'-wave tt-tt scattering data. 

The matrix element for the isospin zero tt-tt scattering 
is expressed as 



A 



/=o 



3A(s, t, u) + A{t, u, s) + A{u, s, t), (14) 



where s,t and u are the Mandelstani variables and the 
amphtude A(s, t, u) is defined in such a way that the 
invariant amphtude for tt^ (pi ) + ttj {p2 ) — ^ tt/j [p^ ) + tt; (p4 ) 
is decomposed as 

5ij5kiA{s,t,u) + 5iu5jiA{t,s,u) + 5ii5jkA{u,t,s) . (15) 

The partial wave amplitude is obtained as 

1 



n/ = 



Ti,=" = |p(s) / dcos6iFL(cos6')^'='J(s,i,u),(16) 



where p{s) = |q(s)|/(167r-/s) with |q(s)| = ^s — 4m2 /2 
and Pl (cos 0) is the Legendre function. 

In the present model, the scattering matrix A[s,t^u) 



is expressed as 
A{s,t,u) 



EsU-^- (17) 



1 

"^Q-niT-m! .,.,,,..„ 

6 ^ — ' ■" ' S — Wis 



It should be noticed that the sum rule (IT^ implies that 



2s 
A{s,t,u) = — , for s<^ml 



(18) 



consistently with the low-energy theorem of the NG 
bosons. After the partial wave and the isospin projection, 
the matrix element for the iso-singlet S'-wave channel is 
obtained as 



Tlz'ois) = Pis) 



g .^TTTTTTTr 



4m" 



E^/.^^ln 



4to: 



nif 



4 



(19) 



r 



where the second term comes from the t- and u-channel 
scalar exchange contributions, and the third term comes 
from the s-channel exchange. In this expression, 7n,r 
stands for the 7r(140) mass. 

Since the rescattering effect should be pro perl y in- 
cluded in the energy region above 800 MeV [ij], the 
above amplitude is applicable only in the low energy re- 
gion. Actually, in the calculation of the Di — > Dtttt decay 
in section |Vl the exchanged energy is y/s < 560 MeV, so 
that, among the exchanged four scalar mesons, the sigma 
meson gives a dominant contribution. One might naively 



eliminate the contributions of fp2, fp3 and /p4 from the 
second and third terms in Eq. (|17p to reduce the num- 
ber of parameters. However, such a truncated amplitude 
cannot reproduce the low-energy theorem in Eq. (|18p ob- 
tained as a consequence of the chiral symmetry. Then, 
instead of naively eliminating the scalar mesons other 
than sigma, we make an expansion of the amplitude in 
terms of s/mj. and m^/mi. (j = 2,3,4). As a result, 
the partial wave amplitude T[z^q{s) is reduced to 



J 



m^ois) = Pis) 



J=2 



Am 



2 i'crTTTT 



In 



4m: 



_J2 ( 5^-1^ -A^(s- Ami' 



-f, 



(™/P^ 



3 4^^ „ 

(m?. )2 



•^SffTTTr" 



(20) 



where dots stand for the higher order terms. 
Making use of Eqs. p^ and (|13l) . we arrive at 



Ti=[;(.) = Pis) 






1 



Ami 



In 



s + TO^ - Ami 



33, 



2 

(TTTTT 



.yCTTTTr 



F^ imlf 



(2s -f Ami 



(21) 



We would like to stress that the amplitude (PT|) depends g^^^^). This low energy reduction is achieved by using 
on only two undetermined parameters Too- and g^TTTr, the relations in Eqs. (1121) and (J13p which are the conse- 
while the one in Eq. (TTOl) includes nine im,f^ , gf^-wTr and quences of the chiral symmetry Thus, for studying the 



Di — >■ Dtttt decay rate, it is enough to determine the 
parameters m^ and gairv 

To study the tt-tt scattering we include the finite width 
effect in the sigma meson propagator. Here, we use the 
modified Breit-Wigner prescription where the width ef- 
fect in the sigma meson propagator is taken to be mo- 
mentum dependent, that is, we make the substitution 



1 



s — m^ s — m'^ + ima-Ta-{s) 



(22) 



The sigma meson width has several expressions in lit- 
erature. In the present analysis, we take the width as 
the imaginary part of the tt-tt loop contribution to the a 
meson self-energy in the linear sigma model 



r(.) 



32tttoo 



35^,^ . /^ Ami 



(23) 



Now, we fit the sigma meson mass ma and the cou- 
pling constant gaTv-rr from the tt-tt scattering data given 
in Refs. [IllllJ] in the low energy region below 560 MeV. 
The best fitted values are obtained as ^ 



ma = 606±9MeV, 
'2 3.48 



\go 



2.16 ± 0.07 GeV, 



X 



d.o.f. 



12 



= 0.29 . 



(24) 



Note that we can determine the absolute value of gcnryr, 
the sign of which becomes relevant in Di — >■ Dtttt decay. 
Here we obtained the ga-irTi in the chiral limit. We expect 
the correction from the inclusion of the pion mass is on 
the order of to^/to^ [J|, which is about 5%. We show 




: [GcVJ 



FIG. 1: Best fitted curve of the real part of the 1 = 
S-wave TT-TT scattering amplitude. The Data are taken from 
Refs. [ll|,|li. 



^ When we use the parameterization of the width of the a 

^>0 I 4-777 rTi 

propagator as T{s) = nS^'" A/l — [J|, the best fitted val- 
ues are m^ = 586 ± 7MeV, |g<T7r7r| = 1.90 ± 0.04 GeV with 
X^/d.o.f. = 7.05/12 = 0.59. On the other hand, when we 

39„,„, /I •."'„ ^g obtain m^ = 597 ± 8 MeV, 

5.23/12 = 0.44. 



: 2.05 ± 0.06 GeV with x^/d.o.f. 



the best fit curve in Fig. [T]with the allowed region of m,^ 
and ga-n-K at Ict level shown in Fig. [51 We conclude that 
our model can reproduce the tt-tt scattering data below 
560 MeV quite well. 




0.605 0.610 

»i,[GeV| 



FIG. 2: Allowed region of nio- and (jcririr at la 



IV. EFFECTIVE LAGRANGIAN FOR 

HEAVY-LIGHT MESONS WITH CHIRAL 

PARTNER STRUCTURE 

In this section, we introduce an effective Lagrangian 
for the heavy-light mesons coupling to the light mesons 
based on the heavy quark symmetry combined with the 
chiral symmetry. For constructing the Lagrangian invari- 
ant under the linearly realized chiral symmetry, we need 
to include the chiral partner to the lowest lying heavy 
quark multiplet oi H ~ {D* , D). In the present analysis, 
we regard the doublet G = {Di, Dq) as the chiral partner 
to the H doublet 0, Q • For constructing the effective La- 
grangian including these fields, it is essential to consider 
the symmetry properties, especially the U{1)a charge, of 
these heavy-light mesons. Since the mass spectra of the 
states in the H and G doublets are consistent with the 
theoretical predictions based on the cq interpretation, it 
is reasonable to regard them as the cq states. Following 
Refs. [i,li,[il[il, we include the H and G doublets into 
the Lagrangian through 

Hl = 4f [G + iH^5] , nR = ^[G- iHj5] .(25) 

In terms of the physical states, the H and G doublets are 
expressed as 



G = 



2 

1 + i! 



HD'^1,75 + D*] , 



(26) 



with v^ being the velocity of the heavy meson. 

Under chiral transformation, the Hl and Hji fields 
transform as 



Hl ^ nLgi, nR-^HRgl 



(27) 



under the axial U{1)a symmetry, Hl and Hr transform struct an effective Lagrangian describing the interactions 
as among the heavy-hght mesons and the hght mesons. In 

the present construction, we only include the minimal 
"Hl — > "Hls '", "Hr — > "Hije*". (28) number of terms which are responsible for our following 

analysis of the Di — >■ Dtttt decay. Then the Lagrangian 
Using the chiral and U{1)a transformation properties of jg -vvrittcn as 
heavy-light meson fields and light meson fields, we con- 



J 



Acavy = ^Tr [HLiiv ■ d)nL] + ^Tr [HRiiv ■ 9)^^] - y Tr [H^Hl + HflH^] 

- ^Tr [M^UlUr + MHrUl] + ^^^Tr [-f^'^M^HLHR - -f''pM'HR%L] , 
4 ZPtt 



(29) 



where A,^,^ and qa are parameters. In terms of the H and G doublets, the effective Lagrangian ([^^ is rewritten as 

•^hcavy 



ixr [-Hiv ■ dH + Giv ■ dG] - — Tr [-HH + GG] 



-^Tr [{M + M^) (GG + HH) ^ i (M ^ A'P) (HG ~ GH) 7^] 
8 

+?:;^Tr [- O'/ - pf^) {GG + HH) 7^ - i {pM + ^Np) {HG + GH)] 



(30) 



r 



V. THE SIGMA MESON STRUCTURE FROM 
Di -^ D-KTT DECAY 

In this section, we study the quark contents of the a 
meson through the Di -^ Dtttt decay. 



Third, we determine the combination g^^ cos 9^^ using 
the Dq — > Dtt decay. The relevant terms in the La- 
grangian (j30l) are given by 



ign 



-)*t 



Cd^d. = ^D<i>D;^ +R.C., 



(34) 



A. Determination of the parameters in the heavy 
meson effective Lagrangian 

In this subsection, we determine the parameters in the 
Lagrangian ((30|) that are necessary for the numerical cal- 
culation of the Di — > Dtttt decay width. 

First, since the kinetic terms of H and G doublets 
have the opposite signs, the A term in Eq. (|30p shifts 
the masses of H doublet and G doublet to the same di- 
rection. Then, in the following analysis, we take A = 
without loss of generality. 

Second, we fix the value of the combination gA cos0^ 
from the partial width for D* -^ Dtt decay. From the 
Lagrangian ((30|). the D*-D-tt interaction is extracted as 



C 



D'D-n 






Ddn<i>D*f''' + H.c. 



(31) 



From this the D* ^> Dtt is expressed as 

{g A cos 9^)^ mjf 



T{D 



* + 



Dtt) 



4^F2 



Ip.I , (32) 



'D' 



where JPttI is three- momentum of outgoing pion and ttih 
is the average mass defined as mn = (Sjtid* -|- mD)/4. 
Using the central value for r£|.+ = 96 KeV [4I we get 



\gA cos 9tj\ — 0.56 



(33) 



From this Lagrangian, the width for Dq -^ Dtt decay is 
expressed as 



T{D* ^ Dtt) 



S mnrnG {gnCos9j,y 
2 2TTmli. 4 



(35) 



where mc is the average mass defined as ma — (3moi + 
"T'D*)/4. Using the central value of the data T{Dq) = 
267 MeV [2| and assuming that D^ decays dominantly 
to Dtt, we obtain 



IffTrCOS^^I 



3.61 



(36) 



With the numerical value of {g^^ cos 9t^\ estimated above 
we predict the Di — > D*tt decay width as a check of the 
validity of the heavy quark expansion applied here. From 
the Lagrangian ([50]) we obtain the relevant Lagrangian 
as 



J^DiD*Tr 



ig^ 



Mt 



D*^D';^ +R.C 



which yields the following expression of the Di 
decay width: 



r(Di -^ D*tt) 



3 TTlHrnG {g-r: COS 9-„Y 

2 27rm|,^ 4 



IPttI 



(37) 

D*TT 

(38) 



Using the central values for the relevant particle masses approach for the H and G doublets, 
we obtain 



r(Di -^ D*tt) = 224.0 MeV. 



(39) 



Comparing this with the data T{Di) = 3831^°^ ± 74 [1| 
we conclude that our prediction, based on the heavy 
quark limit, is consistent with the data. This implies 
that the heavy quark limit taken here is a reasonable 



J 



B. -Di — )■ Dtttt decay 

In this subsection, using the model explained so far, 
we study the Di —> Dirw decay. We show the relevant 
diagrams in the present model in Fig. |3l 






D 



_5 



Di D* 

'WVXW/V/V/v 



D 



Di 



D* 



D 



'WXa r 



(A) (B) (C) 

FIG. 3: Feynman diagrams contributing to Di — !> Dmv decay in our model. 

The matrix element for this process is straightforwardly written as 



Ai = —\/2mGmH 







9iT 2 

: cos I 



-^772 



V'^V ■ p„2 



2^/2 



97r 

; 1 1= ' 

kD*+iTD-/2 2^/2 






V ■ ko* +ir_D./2 



(40) 



where s = (p^j^ + Ptt^)'^ = "*^7r with p^^ and p^j as the momenta of pions in the final states, e^{v) is the polarization 



vector of Di meson satisfying X^po 



by p^. = mD*v^ 



polarization ^M i^ 



'(S/ii/ — Vfj,v^), and the residual momenta ko' and ku* are defined 



k^t and p^, — rriu'V^ + fc^, , respectively. 



We note that, since the maximum value of the tttt invariant mass is restricted as rriT^T^ < mu-^ — mo — 560 MeV, the 
contribution from the lightest scalar meson, the sigma meson is dominant. Along the same method as that was used 
in the tt-tt scattering case, using the sum rules ((T^ and P^ . the matrix element (I40p is reduced as (for a derivation, 
see Appendix [B)) 



M = -V^momH ep(w) < h 



where the mixing parameter h is defined as 



(PtTi +P^2T 



yij-im 2 ; • r t \ 

S — m% + im,aV rj{s) 

g-n COS 6*^ p'i,^ - v^^v ■ p^2 g^ COS Q^ 



( I ^M I ^ hgaTr^\ 






2V2 v-kD-.+ iTo* /2 2 V2 V ■ kDs + iTd* /2 j ' 



COS^^ 



(41) 



(42) 



In the amplitude (HI]), the values of {g-^ cos6'7r| and \gA cos6'7r| as well as those of Too- and l^o-irTrl are phenomenologically 
determined above. Then the amplitude (|4ip depends only on the mixing parameter h. When the physical pion is a 
pure qq state, i.e., cos^tt — 1, then h = 1 implies that a meson is the qq state while h = the pure qqqq state. 

Next, we make the isospin and the partial wave projection with respect to the final two pions and pick up the 
/ = 0, S-wave amplitude. This procedure eliminates the contributions from L = 2, 4, • • • in diagrams (B) and (C) in 



Fig. |3l ^ The isospin projection gives the isospin singlet as 



|7r7r,/ = 0) = --=\TT+{pi)TT-{p2)+TT-{pi)n+{p2)+TT°{pi)TT"{p2)) 

And, the partial wave projection decompose the Di — )• Dttt: decay matrix as 

dr{Di -> Dtttt) .^^^ dVLiDi -> Dtttt) 



drriT, 



E' 



dm.„ 



(43) 



(44) 



where 



^r-(^^-^'^ = 2(2L + l)l 1 



dm-n 



2l64:^T^mD^ \2j Ed 



A Ip^I 



Y, \MLiml^)\\ 



(45) 



polarization 



with IpttI = \/rnu\^ — 4to^/2 as the three momentum of the outgoing pions in the center of mass frame of two pions. 
M^L^Jf^TXTx) stands for the D\ — > ZJtttt decay matrix element with the final two pions in L-wave, explicitly, 



Mhiml^) 



1 /•! 

- / dcos0Pi(cos6')7V{(m7r,r,cos6'), 

2 J ~\ 



(46) 



with A^ given by Eq. (|4ip and 6* being the angle between tt and D. In the following, we only consider the / = 0, 5- wave 
dr(Z?i -^ D{T:T:)j^Q]^^Q)/dmT^^ to study the composition of the a meson. Taking into account of the sign ambiguities 
of g^Trn and g-^cosO-^^ we show, in Fig. 21 the resultant width with h = Q (solid line) and h — 1 (dotted line) for 
four cases: (a) ^o-tttt > and g.^ cosOtt > 0; (b) ga-mr > and g-^ cosOjr < 0; (c) ga-mr < and g^^ cos^jr > 0; and (d) 
gcTTTTT < and g-nCosO-^ < 0. ( Since gyicos^Tr is the overall factor of the amplitude, the sign ambiguity is irrelevant 
to the decay width.) This shows that, for ga-KTr > 0, the width for ft, = 1 is much larger than that for ft, = 0. On the 
other hand, for gaTv-K < 0, the peak position is located around 450 MeV for ft = while above 500 MeV for h — 1. 
This suggests that, when the experimental data become available in the future, the mixing parameter h together with 
the signs of gaTv-K and go-n-Tr could be fitted. 
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VI. A SUMMARY AND DISCUSSIONS 

We study how the composition of the lightest scalar 
meson a affects the Di — )■ Dtttt decay. We construct 
a linear sigma model for the qq and qqqq chiral nonet 
fields, distinguishing two nonets by their U{1)a charge. 
We write down the effective interaction terms among cq- 
type heavy mesons (D, D*, DS(2400), 1)1(2430)) and the 
scalar nonets, in which the existence of the U{1)a sym- 
metry implies that only the qq component of the scalar 
meson couples to the heavy mesons. 

In the light quark sector, we use the most general form 
of the potential for the two chiral nonet fields. Using 
the sum rules among the coupling constants and masses 
of scalar mesons, we express the tttt scattering ampli- 
tude in the low-energy region in terms of only two pa- 
rameters, ttIct and g^TnT- We fit the values of m^ and 



gcnr-TT to the TTTT Scattering data below 560 MeV, which 
is the maximum energy transferred to two pions in the 
Di — > Dtttt decay. We fix the parameters in the effective 
interaction terms to heavy meson phenomenologically. 
We show how the / = 0, 5- wave differential decay rate, 
dr(-Di — > D{tttt) j^Q^L^o) / drnTTTT, depends on the mixing 
structure of the a meson: For ga-Tnr > 0, the width for 
ft = 1 is much larger than that for ft = 0, and the peak 
position moves to the higher energy region for go-u-Tr < 0. 
In the theoretical calculation, we can include a part of 
the final state interaction effect. The final state interac- 
tion among all the three final particles and that between 
one pion and D meson are suppressed by the heavy me- 
son mass therefore their effects are small. For the final 
state interaction between the two pions, in the I — 0, S- 
wave channel, the contribution from the bubble diagrams 
is illustrated as 
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FIG. 4: dr{Di — >■ D{nn)i^o,L=o)/dmTrTv vs Tn-n-Tr with h — (solid line) and h — 1 (dotted line) for gaTnr > and Qt^cosOt: > 
(upper left panel), gcnnv > and g-n cos^^ < (upper right panel), gcr-mr < and g^ cos 6'^ > (lower left panel), gcr-mr < and 
Ptt cos^tt < (lower right panel). 



^ We calculated the D-wave decay width which is less than 5% of 
the S-wave width. 
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The summation of the term.s in the second parenthesis 
in the right hand side of the above equation gives the 
/ = 0, S'-wave S'-matrix for the tt-tt scattering, therefore 
it just contributes a phase factor to the Di — >■ Dnn decay 
matrix element and further, does not change the partial 
width. In this sense, the correction to our results from 
the final state interaction is small. 

In the present analysis, the interaction terms among 
the heavy-light mesons and light mesons is constructed 
in the heavy quark limit with the linear realization of 



the chiral symmetry, therefore the DiDa and D*D7r 
coupling constants are related by the chiral symmetry. 
Since the chiral symmetry is dynamically broken, there 
exists a difference between the DiDa and D*Dtt cou- 
pling constants. To show this chiral symmetry break- 
ing effect, we typically rescale the DiDa coupling con- 
stant by a factor two and illustrate out results in Fig. [5] 
This shows that the tendency of the differential width 
dT{Di — > £'(7r7r)/=o,L=o)/'^"^7r7r does not change. 
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FIG. 5: dV{D\ — !> D{TTn) i=q^l=o) / drnTriv vs m^^ with h — (solid line) and h — 1 (dotted line) for gcnviv > and g^ cos 6.^ > 
(upper left panel), gcrmv > and (;.^cos6.^ < (upper right panel), gcr-mr < and (;^cosS^ > (lower left panel), gc-wTr < and 
g-n cos&TT < (lower right panel). 
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Appendix A: Derivation of the coupling constants 

In this appendix, we derive the relations ^ and (|13p from the general Lagrangian (|4]). Following Ref. [lOj, we write 
the infinitesimal transformation matrices for vector and axial- vector transformations as Ey and Ea, respectively. 
These infinitesimal matrices satisfy the following relations: 



Ey — -Ey, 



e\ = -Ea, TrEA = 3m. 



Then, under the vector and axial-vector transformations, the relevant meson fields transform as 

5v$ = [Ev,<^], 5vS=[Ev,Sl 

5a^ - -i{Ea,S}, 5aS^i{Ea,^, 

5v^' = [Ev,^'], 5vS' ^[Ev,S'], 

5a<^' = -i{EA,S'} + 2iS'T:TEA, 

5aS' = i{EA,^'}~2i^'TvEA. 



(Al) 



Since the effective potential term Vq is invariant under the chiral 811(3)^x811(3)/; transformation implying the in- 
variance under the vector transformation, then we have 



5vVo = Tr 



(9$ dS 



[(S',$)->(5',$')] =0. 



(A2) 



The U(l)yi symmetry is explicitly broken by the anomaly in T^^, so that the axial transformation is expressed as 



SaVo = Tr 



9$ dS 



[(5,$)^(5',<i>')]=<5K,, 



(A3) 



where SVr, arises from the chiral anomaly. Then, using Eqs. ljAip . (|A2|) and ()A3|) as well as the arbitrariness of the 
variations Ey and Ea yields the following generating equations: 



dVo 
dS 



S, 



dVo 



(5,<i>)^(5',$') = o, 



$,S5Ul5,|^} + (5,<i>)->(5',$')-l 



2Tr 



dS' 9$' 



+ 



3ia 



(A4) 



In the following analysis, we use the stationary conditions for the potential Vq given by 

0. 



(S)-. 



dVo 
dS' 



Furthermore, we work in the chiral limit, therefore VEVs of S and S' are proportional to the unit matrix: 
Differentiating Eq. (jA4l) with respect to $, we obtain 



2W2 






2W4 



Differentiation of Eq. (jA4p with respect to $' gives 



2V2 



9$?a$i 



2vi 





-sf 


2..>;( 

m=l 


' d'Vo \ 






to $' gives 










d^V^ \ _ 


= s^ 


3 

m=l 


1 dVo \ 


1 





(A5) 



(A6) 



Along the same method, differentiating Eq. ()A4p twice with respect to S*, S", $ or $' and using 

a/97r" = (a/a$} - a/a$^)/V2, 
a/97r'o = (a/a$f - a/9$'2")/y2. 
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did fa = {d/dS\ + d/dSl)/V2, 

d/dfb = d/dsl 

d/dfc = {d/ds[' + d/ds'i)/V2, 

d/du = d/ds'l 

we obtain the relations between three-point and two-point couphngs in the isospin hmit as 



(A7) 



V2 



d^Vo 



V2 



V2 



V2 



V2 



V2 



V2 



V2 



.dinofdfa/ 
' d^Vo \ 

XdnOdn'Odfa 

^dirOdTT'f^dfb 

' d^Vo \ 
,d{n"fdfj 

' d'Vo \ 
.dinofdfd/ 

dnOdn'Odfc 

dTT^dTT'^dfd 



Vi 



d^Vo 



I a^vb 



dn^dn'Odfa I V2 \ 9(/a)2 / 72 \ 5(7r0) 



i_m^ 



V4 



I 93^0 \ J_/_9^Vb 



\d^^d^'^dhj V2\dfadh 



Vi 



d^Vo 



/^n_ 



IJ^ 



U4 



d{Ti'^fdfal V2\dfadfc/ V2\a7rOa7r'" 



V4 



-Vi 



\d{7r'0 fdfb/ V2\dhdfc 

dT:^dT:'^dfJ V2\dfadfc/ ^\dT:^d^'^ 

d^^dTT'^dfdl V2\dfadfd 



Vi 



■ v^ 



d^V^ 



^din'oydfci 
,d{7T'0fdf. 



_ 1 / a^^o \ 1 / d^Vo 



V2\d{fc)y V2\d{TT'oy 

1 / d^Vo 
V2\dfcdfd 



(A8) 



And, differentiating Eq. (jA4p three times with respect to S,S',^ or $', the relations between four-point and 
three-point coupling in the isospin limit are obtained as 

d^Vo \ , I aVo 



V2 



V2 



V2 



9(7r0)4 
/ d^Vo 



Vi 



_ , 1 / d^Vo 
9(7rO)3a7r'V ^/2 \a(7r0)29/„/' 



Vi 



5Vo 



V2 



\dTT°d{Tr'°fdTTO ^ 

d^Vo 



V4 



d{TT'°ydTT° 



Vi 



dTT'Odin^ydn'O / 

d^Vo \ _ 1 / g'^o 
dn^din'^f / V2 \d{7r"^ fdfa 



1 / d^Vo \ 1 / d^Vo 

V2\d7T'0d7r0dfa/ V2\d{7T°fdf, 



g 1 / d^V, 
V2\dTr0dTr'0df, 



, 1 / d^Vo 
'V2\d{n'°fdfc 



(A9) 



The above relations are re-expressed in terms of the mixing angles for the scalar and pseudoscalar mesons in Eqs. ^ 
and (O , the pion decay constant and the VEVs of M and M' fields as 



V2^^ 



Fjr cos Ott + Ftt sin O-n 



7^Y.^F^)'^^U-^) 



aj, 



J = l 



V,^- 



-Fn sin 6'7r + -Ftt cos 0^ 



;tE(^-)^(^^"')^^' 



(AlO) 



i=i 



where {Ft^)i = F^, (F^)2 = Ft^ are the decay constants of 7r(140) and 7r(1300), respectively, and the mixing matrix [/^ 
is defined by 



U^ = 



(C^7r)la (t^7r)l6 
{UTr)2a (t^7r)26 



cos( 
sin 6 



- sin Ot, 

COS0, 
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Note that, in the chiral Hmit, only the NG boson can couple to the axial vector current, so that F^^ = 0. Then, writing 
the mass matrices as 



/ d^v \ 



d^v 



where subscripts A and B run from a to 6 and C and D from a to d, we have 



1 / d^V \ 

1 '' 



V2 



d'Vn 



A— a 
b 



di7T")ad{7rO)Adfb 



^F^ 5^^-)^^ \a(7r")^9(7r"),a/„ 



{rn})ab, 



{mf)ac - {m^)ab 



1 / d"^V \ 



V2^" t''^'""^\9i^")ad{7r°)Ad.fc 



1 '' 

1 '' 



d'Vn 



A=a 
b 



53^0 



9(7r"),a(^o)Aa/d 



(-?.).c - (m^).., 



{rh})ad, 



A=a 
b 



a^Fo 



V2^''S^^"^'^\a(vr"M(7rOM/. 



■y = {m})cc- {ml)bb, 



1 / 9 y \ 

~^^-^ y'(^7r)lA ( a/ n\ a^ n\ o i- / ^ ("^f)cd 

\/2 f^ \din")hdiTrO)Adfd/ ^ 



Summing these equations with mixing matrices t/jr and Uf we get the following physical coupling 



1 



--F^ 



d^Vo 



V2 ^ \d{7r0),d{7r0),df, / y/2 



-jn^- E E (u.)iA(u.)kB(Uf),c 



d^V. 



A,B=a,b C=a,b,c,d 
d 
^(Uf)jC [iU^)ka{m'j)aC + (U^)kb{m'j)cC 



di7r")Bd{7r")Adfc 



C=a 



Xl(^'')fe-4 [iUf)jaiml)aA + iUf)jciml)Ab] , 



A=a 



where we write the physical pseudoscalars as (7rp)i = tt^ and {Trp)2 = 'n'p ■ 
So that, in the case that the pseudoscalars are neutral, we have 



9f,7 



d^V 



a«)i9«)i9/, 



b b d 



EEE(^-)i-4(C^-)ib(^/) 



JC 



A— a B—a C—a 



d^'Vo 



d{7r0)AdiTT0)Bdfc 



As a result 



9f,7:7T 



J2 E iu^)Mu.)iB{Uf) 



JC 



A,B=a,h C=a,b,c,d 
2 / ^ n2 



d^Vo 



di7r")Bd{7r")Adfc 



— (ml) [{U.)la{Uf),a + {U.)lb{U}),c 



(All) 



(A12) 



(A13) 



(A14) 
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Here we have considered that in the chiral Umit, the hght pions are massless. Therefore, from this equation, we 
obtain the following relations 

4 2 4 



H) 



2 A^ p2 



V2 






-1 "^f: 



Similarly, in the the chiral limit, we can derive that 4-7r^ coupling constant satisfies 



Q-K-n 



d^V 



a«)i9«)i9«)i9(<)i/ ^2^ 



Y^ 2 



H 2 



(f/0la(f//),a + (C/.)l6(f//),. 



which yields 



by using Eqs. (|ATl [Aiell . 



"" 0^ 1 

,^2 — ^BiTTT-mn 



-1 "^'h 



(A15) 



(A16) 



(A17) 



Appendix B: Reduction of the matrix element (|40|) 

Here, we provide the reduction of the heavy meson decay matrix element (|40p . Its first term can be rewritten as 









i=2 



rrir + inif.Tf. (s) 



Considering that the momentum carried by the exchanged scalar meson is smaller than the heavier scalar meson mass, 
i.e., s < mi, ,{i — 2, 3, 4), we make the expansion 






i=2 



mi, + im,f.rf.{s) 






i=2 



1 



mj^ -^rrifTf^is) 



And also concerning that Tj. <^ "m-f-, (i — 2, 3, 4), we can simplify the above relation as 



E^/.' 



(PtVi +P^2y 



^^2 s-m2^ +im/T/^(s) 



With respect to the relation 



fTr-l\ _ \^ „ (P^l +P^2Y ITT-I 



—^{Uf )« = T^cos^,, 



1 



2 ' I 2 
rrif, \ rrir. 



— rrif 

=1 fi 



F^ 



and neglecting the terms of and higher than 0{s/m,i,), {i = 2,3,4) we obtain 

E{PtT1 ~^ PtT2) ITT-1\ ( I ^M I V2 „ gaiT-K ,jj-l 



i=2 



s — mi, + imf. i- fi(s) ■' 



mi: 



iuy%i 



(Bl) 



(B2) 
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Then, we finally arrive at the sum rule for the scalar meson exchanging contribution 

4 



(Ptti +P7T2T 



E KPtvi -^ Pti2 r nT-i\ 



(PtTi +P^2)^ 



TTTT O 



-(C/7^)al 






(B3) 



It should be noticed that the second term of this sum rule arising from the heavier resonance contribution and only 
the first term leaves once the heavier resonances are neglected from the beginning. 
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